Modulated vortex states in Rashba non-centrosymmetric superconductors 
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Vortex lattice structures to occur in Rashba non-centrosymmetric superconductors under a mag- 
netic field parallel to the basal plane are studied by assuming the s-wave pairing and taking account 
of both the paramagnetic and orbital depairings. A high field vortex lattice with a modulation 
perpendicular to the field and of Fulde-Ferrell-Larkin-Ovchinnikov type, occurring in the limit of 
vanishing £ (spin-orbit coupling normalized by Fermi energy), is found to be replaced with increas- 
ing £ by another new vortex lattice with a modulation of the gap amplitude and induced by the 
absence of inversion symmetry. A correlation between a structural transition and an inflection of 
the _ff C 2(T')-curve is discussed in relation to H C 2 data of CeRhSi3. 

PACS numbers: 



Recently, superconductivity in materials with no in- 
version symmetry has been a focus of intense interest 
[3, 0] . A spin-orbit coupling due to the absence of an in- 
version symmetry may lead to coexistence of spin-singlet 
and triplet pairing states and result in newphenomena 
stemming from the pairing symmetries 0, [4[ . Another 
interest is a consequence of the non-centrosymmetry in 
superconductivity in nonzero magnetic fields. An inter- 
play of the Zeeman effect and a splitting of Fermi surface 
(FS) due to the spin-orbit coupling was shown to, in fields 
(H _L c) parallel to the basal (a-b) plane, lead to a he- 
lical phase modulation Q of the superconducting order 
parameter A(r) similar to the Fulde-Ferrell state [|| in 
a Pauli-limitcd superconductor as far as the densities of 
states Nj on the splitted two FSs, denoted by j = 1, 2, 
are different from each other. Later, the mixed state was 
studied by invoking the Pauli limit [jj with no vortices, 
and, just as in the centrosymmetric case, the presence of 
a modulation of the amplitude |A| of Larkin-Ovchinnikov 
(LO) type [8[ was argued which is separated by a single 
transition from the phase-modulated state in lower fields. 
However, it is natural to expect appearance of a different 
superconducting phase from those in the centrosymmet- 
ric case. 

In this work, mixed or vortex states occurring in 
Rashba non-centrosymmetric materials jl| in H 1 c are 
investigated in clean limit by assuming the s-wave pairing 
and incorporating both the paramagnetic and orbital de- 
pairings. Based on existing estimates |4| of the spin-orbit 
coupling energy £ in real non-centrosymmetric materials, 
the relation Ep f ^> T c between relevant energy scales 
will be reasonably assumed, where Ep is the Fermi en- 
ergy, and T c is the zero field (H = 0) superconducting 
transition temperature. Then, the key parameter mea- 
suring the broken inversion symmetry in superconductiv- 
ity is SN = \Ni - N 2 \/(N 1 + N 2 ) which is of the order 
of £ = Q/Ep rather than k^Tc/Q. The main focus in 
H _L c, where the paramagnetic depairing is effective, is 
to understand effects of the absence of inversion symme- 
try on a modulated state of LO type in higher H and 
at lower temperatures. Our findings are summarized as 
follows: In the specific limit, SN — > 0, the LO-like state 



with a modulation perpendicular to H occurs as a result 
of the paramagnetic depairing and is transformed with 
decreasing H into an Abrikosov lattice through two first 
order structural transitions (FOSTs). However, changes 
of the lattice structure obeying the shift of FSs induced 
by the Zeeman energy do not occur. Once a small but 
nonzero SN is taken into account, however, the LO-like 
state is pushed up to higher H, and a new intermediate 
phase appears accompanied by an amplitude modulation 
parallel to the direction of the helical phase modulation 
as a consequence of the lack of inversion symmetry. As 
\SN\ is increased so that the LO-like state disappears, the 
high field side of the H-T phase diagram is entirely occu- 
pied by this new modulated state. By comparing these 
results with those in the Pauli limit [3] , essential roles of 
the orbital depairing in this issue are emphasized. 
We start from a BCS-like Hamiltonian 

Ubcs = X! c k j;il ( £ k<5 Sl , S2 + (Cgk + MsH) • <t S1}S2 ) c k , S5 

k,si,s 2 



(1) 



where fJ-sH is the Zeeman energy, £k denotes the bare 
band energy, and gk is the unit vector (k X I) with 
z perpendicular to the basal (a-b or x-y) plane. A 
Ginzburg-Landau (GL) free energy F will be derived 
within the mean field approximation and in type II 
limit with no variation of flux density. As performed 
elsewhere [![, diagonalization of the quadratic term 
leads to a splitting of FS into two branches (j = 1, 
2), and, consistently with this, the pair- field operator 
* q = -iEkKW-k + q/2(3 c k + q /2 a may be rewrit- 
ten in terms of the field operator a^j on each band as 
EkEj ; /j( k ) a -k + q/2j -ak + q /2j 9], where fj(k) satis- 
fies \fj\ = 1, and its detailed expression is not neces- 
sary within the weak-coupling BCS theory. In reaching 
the latter expression, a small contribution due to inter- 
band interactions of O(e^) to F was neglected, where 
Eh = Max(/is_ff, T c )/£. By replacing the interaction 
term in eq.(l) by J d 3 r(|A| _1 1 A| 2 - & A - *A*) in terms 
of the order parameter A, expanding the free energy in 
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A results in F . Its quadratic term is given by 



F 2 = / tfVA*(r) 



|A| 



dp 



27rTN(0(p;U Q )) 
sinh(27rpT) 



A(r), 



where N = Ni + N 2 , ( ) denotes the average over the 
momentum k on FSs, and 



o(p-,n Q ) = i£ 



cxp(-ipcrv x Q)(cos(2p B Hk x p) (2) 
cxp(-icrpv • IIq) 



a=±l 

+ la 5N sm(2p B Hk x p)) 



under H parallel to y, where IIq = — iV + 2eA — Qx, 
v = v_l + v z z is the velocity vector on FSs, and a possi- 
ble modulation wave vector Q = Qx was introduced (see 
below). The ^-dependence accompanying p B H implies 
that the paramagnetic depairing does not work effectively 
in the free energy terms associated with the gradient par- 
allel to H. Thus, in contrast to the centrosymmetric case 
[ToL [Tlj . A(r) solutions modulating along H never ap- 
pear in Rashba superconductors. In H <C izp(O), where 
H P (0) = A(0)/V2MB is the Pauli-limiting field, the fac- 
tor prior to exp(— iupv • IIq) is approximated by unity 
if Q = Qo = 2p, B H5N/\v 1 _\ [|, and then, A is diag- 
onalized in terms of Landau levels (LLs) with even in- 
dices, A = J2 n °2» <f>2n(r - Qor 2 H z\0), where r H 2 = 2eH, 
and 4>o (r 1 0) is an Abrikosov lattice solution in the low- 
est (n = 0) LL and is one of magnetic Bloch states [Icj . 
(f> n (r\r ), in the same gauge A = Hzx. Note that, us- 
ing the property «^„(r|r ) = exp(iz x/r 2 H ) </>„(r + r |0) 



!or 2 H z). 



helical state [5j, A = 
Therefore, the heli- 



[IOJ, we have the form of the 
exp(iQ x) S n a2n^2n(r| 

cal phase modulation [5[ in H <C Hp merely corresponds 
to a uniform displacement of the vortices. By contrast, 
in H > Hp(0), the 5N dependent imaginary terms can- 
not be erased from the factor prior to exp(— iapv -Hq) in 
eq.(2) any longer and result in a A-solution represented 
by both of even and odd LLs. This hybridization of even 
and odd LLs in equilibrium is a new feature stemming 
from the lack of inversion symmetry and, as is shown 
below, leads to a new vortex state. 

Hereafter, we use a bare FS of the form 10] of a 
corrugated cylinder extending along z and with a co- 
herence length anisotropy 7 = £, x /£,z = 6, where £ M is 
the coherence length in ^-direction. However, no layer 
structure is assumed in real space. Quantitatively sim- 
ilar results should follow in terms of an ellipsoidal FS 
with a similar 7-value as far as H _L z. The paramag- 
netic depairing is measured by the Maki parameter in 



* (2) 



7.lp B H^>(0)/(2nT c ), where H^'(0) 



(2D) 



IS 



H|| 

the 2D orbital-limiting field. Note that the argument of 
the trigonometric functions in eq.(2) is proportional to 



a$h, where h = H/H^'(0) 
high field state in low T is, as in various centrosymmetric 
systems [III [l2j], described by higher LLs. 



2D), 



At least in > 2, 



the 



In examining the details of H-T phase diagrams by 
including higher LLs, the best approach will be to diag- 
onalizc F2 in terms of A = ^ n a n </>„(r — Qr 2 H z\0) with 
the consistent minimization in Q, dF/dQ — [9( im- 
plying that the total current is zero, where </>o( r |0) = 
J2 p exp[ikpx - 7(2 + kp r 2 H ) 2 / (2r 2 H ) + i(k 2 r 2 H p 2 cot0/2) ], 
and the parameters k and 9 determine the vortex lattice 
structure. Due to inclusion of numerous LLs, however, it 
is not tractable to perform the Q-minimization exactly. 
As an alternative approach, we first choose Q — Qo and 
try to cure the deviation Q — Qq by incorporating more 
numerous LLs. In this procedure, a finite uniform cur- 
rent due to the deviation Q — Qo should be cancelled by 
the corresponding one [l(| arising from the even-odd LL 
mixings [13| . Besides, inclusion of more LLs are needed, 
even when SN = 0, with increasing [Tll.[l2|. Thus, a 
truncation of LLs might lead to an ambiguity in numeri- 
cal results with increasing \ SN\ or aj^ . However, we have 

verified in the range a M < 5 and for < SN < 0.2 that 
vortex lattice structures following from eight LLs (n < 7) 
are essentially unaffected by adding higher two LLs fur- 
ther. Below, we will focus on A expressed by the eight 
LLs with Q — Qo. Then, by diagonalizing F2, the result- 
ing superposition of LLs with the lowest energy is chosen 
at each (T,H) as a candidate of A in equilibrium. We 
note that the details of FSs are encoded as the relative 
weight between LLs in A. 

Next, after substituting the resulting candidate of A 
into the quartic term F4 of F, the parameters (k,8) in 
0o(r|O) are determined by minimizing F& at each (T,H). 
However, a full evaluation of F4 is a hard task because 
of spatially nonlocal corrections [ll| in F4 and of the 
inclusion of numerous LLs. Fortunately, since roles of 
the nonlocal corrections tend to be reduced with increas- 
ing [13], it is proper in the present case to replace 
-F4 derived microscopically by a spatially local conven- 
tional one, F^ioc = C4J G? 3 r|A(r)| 4 , where the constant 
C4 is derived microscopically. Further, we have verified 
that F4 > by assuming A expressed by a single LL 
or by a couple of LLs as the candidate, suggesting that 
the i? C 2-transition is of second order in Rashba supercon- 
ductors. For these reasons, we have examined the lattice 
structures using i^joc with a positive C4. Note that, to 
determine the stable structure at each (T,iJ), the H and 
T dependences of C4 are not needed. 

Let us start with explaining the resulting phase di- 
agram in the limiting SN — case where Q = @ 
(Fig.l). We find that, through diagonalizing F 2 , H c2 (T) 



FIG. 1: (Color online) A h-t phase diagram for q^' = 2.85 
and SN — obtained by including LLs with < n < 7 in 
A, where t = T/T c o- Each red (blue) solid curve denotes an 
FOST (H C 2(T)). Insets show real space patterns of |A(i,«)| 2 
at the indicated (t,h) points. The axis of cylindrical FS is par- 
allel to z, while Q || x. The Pauli-limiting field corresponds 
to h = 0.496. 
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FIG. 2: (Color online) Corresponding ones to Fig.l for (a) 
SN = 0.003 and (b) 0.2. The state below FOST2 is the tri- 
angular lattice everywhere. 



FIG. 3: (Color online) A h-t phase diagram for \SN\ = 0.003 
obtained in terms only of the n = and 1 LLs of A and 
a<$ = 2.85. 



for SN = is separated by an FOST (FOST1) into one 
branch above FOST1 determined by only odd LLs and 
another one below it due only to even LLs. We have ver- 
ified that this division into the odd and even branches 
for SN = remains valid at least in the range < 7. 
Above FOST1, the state below H C 2 is the LO-like phase 
formed by only odd LLs, while it has no additional nodal 
lines below FOST1 where only even LLs participate in. 
As the insets (A) and (B) show, the nodal lines in the 
LO-like state are removed in passing through FOST1 by 
decreasing H. In any state we obtain, the periodicity in 
z of |A| appears as a periodicity of vortex positions. The 
energy of (A) is quite close to but slightly lower than that 
of a similar structure with nodal lines perpendicular to 
Q. Namely, the orientation of nodal lines is controlled by 
interactions between the vortices and unaffected by the 
details of FSs. In fact, the spacing between neighbor- 
ing nodal lines is not of 0(\v±\/(1bH) but of O(rjj), 
reflecting the details in real space rather than in the k- 
space. The lattice (B) is continuously transformed into 
a deformed square lattice (C) with a smooth rotation of 
the wreckage of the LO nodal lines. Figure 1 includes an- 
other weak FOST (FOST2) in lower H: The lattice (C) 
dis ■continuously changes through FOST2 into the famil- 
iar triangular one (D) with decreasing H. Since the odd 
LLs do not appear around FOST2, the origin of FOST2 is 
the mixing of even higher LLs with n < 6 becoming more 

(z) 

important for larger a\, . This transition on the lattice 
symmetry might correspond to that suggested in Ref. [f| • 
However, it was not considered there [a, @] consistently 
with the appearance of the LO-like state in higher H in 
the same phase diagram. Besides, the fact that some 
lattice structures in Fig.l have no reflection symmetry 
in z suggests a possibility of two lattice domains in this 
artificial case with SN = 0. 

Figure 2 shows examples of the H-T phase diagrams 
in SN case where the even and odd LLs are mixed 
with one another. Even a small enough \SN\ (< 0.003) 
divides the intermediate phase in Fig.l surrounded by 
FOST1, FOST2, and H c2 (T) into two regions separated 
by another FOST (FOST3): Once SN becomes nonzero, 
a new striped phase, (B) in Fig. 2 (a), begins to appear 
near the three transition lines surrounding (B) and (C) 
of Fig.l. Due to a slight increase of \SN\, the growing 
new phase pushes the region including (D), similar to 
(C) of Fig.l, down to lower temperatures (see the arrows 
in Fig. 2 (a)), and the nodal lines of the LO-like phase 
(A) disappear reflecting the even-odd LL mixing. The 
new striped state consists of rows of vortices, separated 
by two neighboring stripes perpendicular to Q on which 
|A| is nonvanishing, and the inter stripe spacing is the 
structural period parallel to Q. With decreasing h, this 
new phase shows a continuous crossover to a square lat- 



tice, (C) in (a), just above FOST2. Besides, a further 
increase of \SN\ not only expands the region of the new 
phase and the square lattice to shift FOST2 to lower h 
but also elevates FOST1 up to H c2 (T) so that the LO- 
like phase is lost. Since \SN\ oc C, the expansion at larger 
\SN\ of the new phase supported by the nonzero Q || x is 
reasonable. Through a comparison with Fig.l, it is seen 
that an increase of \5N\ enhances H C 2 and compresses the 
structure along z. Besides, Figs. 2 (a) and (b) also show 
that, irrespective of\5N\, FOST2 merges with an inflec- 
tive point of H C 2{T) in each figure at which d 2 H C 2/dT 2 
changes in sign. This feature may become relevant in 
experimentally verifying vortex lattice structures found 
here. 

To clarify implication of the above result further, it 
will be compared with Fig. 3, in which only the two low- 
est LLs (n = and 1) are kept, and A is expressed by 
Aoi = a (f) + a\4>x- When SN — 0, the situation is es- 
sentially the same as that in a centrosymmetric system 
near H C 2 ll| : The two LLs do not coexist in equilibrium, 



and a single FOST separating two states with different 
structural symmetries occurs. The low field state is the 
triangular vortex lattice with a\ = 0, while the high field 
phase expressed only by the n — 1 LL is an LO-like phase 
composed of an alternation of a vortex row and a nodal 
line [12| . The situation is similar to the behavior occur- 
ing through FOST1 in Figs.l and 2 (a), where the state 
(A), similar to this n = 1 LL state, is dominated by the 
n = 1 LL. However, it is difficult to understand the global 
phase diagram by concentrating on Aoi. In fact, as al- 
ready noted, the square to triangular structure transition 
at FOST2 in Fig.l is not found without including even 
higher LLs (n > 2), and Fig. 3 shows that the new striped 
phase and FOST2 are not obtained within Aoi. 

The fact that the LO-like state with periodic nodal 
lines is pushed up to the high H and low T corner of the 
H-T phase diagram to disappear with increasing |5iV| is 
in conflict with the corresponding result [7J obtained with 
no vortices in which the LO-like phase tends to disappear 
away from the H C 2 line. In fact, the results in Fig. 2 sug- 
gest that diappearance of the LO-like state is facilitated 
by the growth of the new vortex phase, absent in the 
Pauli limit Q , with an amplitude modulation parallel to 
Q, implying that inclusion of vortices is indispensable in 
the present issue. On the other hand, the H C 2{T) curve 
is similar to the result in the Pauli limit except the fact 
that it is not divergent but flattened in T — ► limit. 

Although we have neglected a p-wave pairing compo- 
nent [151 ] to be induced by the s-wave pairing due to the 
nonzero £, the present results are still valid if the p-wave 
pairing interaction is repulsive. The same thing holds in 
the case dominated by the p-wave. In the case, presum- 
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ably applicable to CePtaSi, where both the singlet and 
triplet pairing channels are attractive and equally impor- 
tant, however, another even-odd LL hybridization due to 
a field-induced coupling between these two pairing states 
dominates over the hybridization appeared in eq.(2) (l6j 
and may change the pictures mentioned above. Results 
of other extensions will be reported elsewhere. 

In conclusion, the vortex lattice structures and H-T 
phase diagrams to occur in Rashba non-centrosymmetric 
superconductors in H _L c have been examined, and even 



a small value of normalized spin-orbit coupling £ is found 
to destroy the LO-like state and replace it by a new 
striped vortex lattice. The obtained results, including 
the square to triangular lattice transition close to the 
inflective point of H C 2(T), are visible in materials with 
standard values of the Maki parameter and anisotropy. 
The present results are relevant to CeRhSi3 in H 1 c 
where the inflective iJ C 2(T)-curves have been found p7| . 

We are grateful to S. Fujimoto and Y. Yanase for dis- 
cussions and to T. Saiki for supplementary calculations. 
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